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We propose a practical way of circumventing the sign problem in lattice QCD simulations with 
a theta-vacuum term. This method is the reweighting method for the QCD Lagrangian after the 
SI/a (3) 8) Ua(1) transformation. In the Lagrangian, the P-odd mass term as a cause of the sign 
problem is minimized. Additionally, we investigate theta-vacuum effects on the QCD phase dia- 
gram for the realistic 2+1 flavor system, using the three-flavor Polyakov-extended Nambu-Jona- 
Lasinio (PNJL) model and the entanglement PNJL model as an extension of the PNJL model. The 
theta-vacuum effects make the chiral transition sharper. We finally investigate theta dependence 
of the transition temperature and compare with the result of the pure gauge lattice simulation with 
imaginary theta parameter. 
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1. Introduction 

The existence of instanton solution requires the QCD Lagrangian with the theta vacuum: 

se = ZqArv D v+ m f)if + - id ^2 £ ^ P W P , (i-D 

in Euclidean spacetime. Though the angle can take any arbitrary value theoretically, experimental 
measurements of neutron dipole moment give the upper limit of theta, |0| < 1 9 [[TJ, g]. Why 
should be so small? This long-standing puzzle is called the strong CP problem. 

Since the upper limit is determined only at zero temperature, the behavior is nontrivial for finite 
temperature. Hence the first-principle lattice simulation is needed, but it has the sign problem for 
finite 0, since the action is complex there. For this reason the lattice simulation is only performed 
by Taylor expansion around = or analytic continuation from the imaginary region where the 
action is real j3|]. On the other hand, after making £/a(1) transformation 

q = e 7i ^q\ (1.2) 
dependence appears only through the mass term 

m f (d) =rnfCos(d/Nf) + rn f iy 5 sm(d/Nf), (1.3) 
in the transformed Lagrangian 

Sf = £q f (YvD v + m f (0))q f + t^ v F* v . (1.4) 
/ * 8 

The P-odd mass term including iy=, makes the fermion determinant complex. 

We propose the following approach in order to circumvent this sign problem^]. Performing 
SUa(3) 8) Ua(1) transformation, 

q u = e iys ^ql , q d = e iy ^q' d , q s = qj, (1.5) 

one can find that Lagrangian ( |l.ip becomes the following form 

JSf= £ q/^mqZ + qJ^qJ + ^F^F^ (1.6) 

l=u,d q S 

Jtl{Q) = y v D v + m/cos(e/2)+m//y 5 sin(e/2), (1.7) 
J£ s = Y v D v + m s . (1.8) 

Here 6 dependence appears only in the light-quark-mass term and the sign problem is induced by 
the P-odd term. However, the scale of this P-odd term is much smaller than Aqcd and hence this 
term is expected to be negligible. If the P-odd mass term is neglected in the reference theory of the 
reweighting method, the expectation value of operator is obtained as 

= j ®A 0'(det^,'(d)) 2 det^ s e- Ss , (1.9) 
'(d) = YvDv+mi cos (0/2), (1.10) 



(deL#/(0)) 2 ' 



(1.11) 
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with the gluon part S g of the QCD action. If this reference system is good, we can expect 6 ' « & 
and calculate (G) with good accuracy. 

In order to justify our proposal, we use effective models and investigate vacuum condensates. 
First, we compare two results with and without the P-odd mass term and examine the effect of the 
neglect. Secondary, we investigate the phase structure in the Q-T-\i space and 6 dependence of 
transition temperatures. This work is mainly based on the Ref. [|J]. 

2. Model setting 

We use the three-flavor Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model 

8 

££ = q f {y v D v + m f -y A }x)q f -G s Y, [{qKqf + (qiJsKq) 2 } 

a=0 

+G D [detq f (l-Y 5 )q f + detq f (l + r 5 )q f ]+^(TMA},^*[A]), (2.1) 

where D v = d v — i8 v ^A^ and the Gell-Mann matrices X a act on the flavor space. The three-flavor 
quark fields q = (q u ,qd,<Js) have masses riif = diag(m„,m^,m s ), and the chemical potential matrix 
jX = diag(jU,jU,ju) is denned with the quark-number chemical potential fi. Parameters Gs and 
Gd denote coupling constants of the scalar-type four-quark and the Kobayashi-Maskawa-'t Hooft 
(KMT) determinant interaction respectively, where the determinant runs in the flavor space. 

The KMT determinant interaction breaks the E/a(1) symmetry explicitly. We use the Polyakov 
potential % of Ref. 



W(T,&]A],&*]A]) 



- C ^-<S>*<S> + b{T)\n (1 - 6<I>*cJ> + 4(<£ 3 + <I>* 3 ) - 3(<J>*0>) 2 ) 



,(2.2) 



( ^)+^(^) 2 , b(T) = bJ^\\ (2.3) 



as a function of the traced Polyakov-loop <J>. 

The four-quark vertex Gs is originated in a one-gluon exchange between quarks and its higher- 
order diagrams hence the Gs can depend on <£. We simply assume the following form for Gs f§], 

G s ->G S {<$>) =G S [l-ai<£*cJ>-a 2 (<£ 3 + <£* 3 )] (2.4) 

which preserves the chiral symmetry, the charge conjugation (C) symmetry and the extended Z3 
symmetry The effective vertex Gs is called the entanglement vertex and the PNJL model with 
this vertex is the EPNJL model. It is expected that dependence of Gs will be determined in future by 



the accurate method such as the exact renormalization group method Q10|, |12| ]. The parameters 
Oli and ai are fitted to (011,012) = (0.25,0.1), to reproduce the result of degenerate three-flavor 
LQCD with imaginary chemical potential [[D]]. 

The EPNJL model has good consistency with lattice results. For the transition temperature, 
the PNJL model is good for the deconfinement transition but overestimates the lattice data for the 
chiral transition. However the EPNJL model well reproduces both of the lattice data. Additionally, 
properties of QCD in the pure imaginary potential region are important. In this region, it has a 
periodicity in Im(ju)/r called the Roberge-Weiss (RW) periodicity 01 . The PNJL and EPNJL 
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models also have the periodicity. At boundaries of the RW period, there is a critical point called 
the RW endpoint. Recently, quark-mass dependence of this endpoint was investigated by lattice 
simulations [15, |T^]. The EPNJL model is successful in reproducing the LQCD result, while the 



PNJL model cannot reproduce this property []13|]. 

6 dependence of the EPNJL model is introduced through the KMT interaction: 



J2f = q f (y v D v +m f -y A fl)q f -G s {<S>) £ [{qKqf + {qiYsKq? 



a=0 



e w detq f (l-Y5)qf + e- ie detq f (l + y 5 )q f , + W(T,<t>\A],®*\A]). (2.5) 



-Gd 



Performing the chiral transformation (|1.5[), one can get 



if = q}(y v D v +m f (6) - y 4 fl)q} - G s (0) £ [(q 'Kq '? + (q 'ijsKq '?] 

+G D [detq}(l - y 5 )q} ! + d e tq}(\ + y 5 )q},] +<&(T,&\A],&\A]), (2.6) 
mi (d) = m,cos(d/2)+miiy 5 sin(d/2), (2.7) 
m s (d) = m s . (2.8) 

Here 6 dependence appears only in the mass term for light quarks. The P-odd term has a much 
smaller scale than Aqcd- 

Since P symmetry is broken at finite 6, we consider P-even and P-odd condensates, 



Of = (?/?/) 

nj = (qfiysqj) 



(2.9) 
(2.10) 



with f = l,s and assume isospin symmetry (a u =0^ = 0/, T] i( = T] f/ = r\ t ). The vacuum conden- 
sates X = aj , rjj /,<£ and are determined by the stationary conditions, 



da 
dx 



0, 



(2.11) 



where Q. is thermodynamic potential calculated with the mean field approximation. 



3. Numerical results 

Figure |l| shows 6 dependence of condensates at T = fx = 0. Since O = at zero temperature, 
the PNJL and EPNJL models give the same result as each other. In Fig. [j](a), the P-odd mass 
is taken into account. Solid and dashed lines show the P-even condensates for light and strange 
quarks, while dotted and dot-dashed lines correspond to the P-odd condensates. Since P-odd con- 
densates are much smaller than P-even condensates, P-odd mass is expected to be negligible for 
this case. In Fig. [T](b), the P-odd mass term is neglected. P-odd condensates (T]/', ?]/) become zero, 
but P-even condensates (a/', O s ') are not affected by the neglect. Therefore it is expected that vacua 



with and without the P-odd mass term are similar to each other and the reweighting method (1.9) 
works well. 

Figure show phase diagrams in the T-pL plane as a function of 6 obtained by (a)the PNJL 
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Figure 1: dependence of the order parameters at T = jJ. = in the EPNJL model. Panel (a) shows a result 
with the P-odd mass and panel (b) corresponds to a result without the P-odd mass. 




Figure 2: Phase diagram of the chiral transition in the [i-9-T space. Panel (a) shows a result of the PNJL 
model and panel (b) corresponds to a result of the EPNJL model. 



and (b)the EPNJL model. Solid lines show the first order chiral transition and dashed lines corre- 
spond to the chiral crossover. Hence, point A is a critical endpoint at 6 = 0. The endpoint slightly 
moves to smaller /I and higher T when 6 is increased from to n. In the PNJL model, the critical 
endpoint does not disappear even for 6 = %. However, for the EPNJL model, this movement is 
much faster and there is no critical endpoint at large 6. This means a possibility that the cosmic 
evolution is changed at QCD epoch by the first order transition if 6 is large. 

Figure |3| shows 6 dependence of transition temperatures at /I = 0. Dashed and dotted lines 
show EPNJL model results for the deconfinement and the chiral transition temperature, respec- 
tively. Solid line shows a result of lattice simulations^]: 

||^ = i-/? e e 2 + o(e 4 ), (3.i) 
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Figure 3: 6 dependence of transition temperatures at jj, = 0. The solid line shows the lattice result for 
the deconfinement transition [|[. Dashed (dotted) lines correspond to result of the EPNJL model for the 
deconfinement (chiral) transitions. 



Re = 0.0175(7). 



(3.2) 



The coefficient Rg has been determined by lattice simulations of pure Yang-Mills theory with 
imaginary 6 parameter, and the constant T c (0) is fixed to that of EPNJL model. Compared with 
the lattice result, 6 dependence of model result is much smaller. This result shows that lattice 
simulations with dynamical quarks are crucial for theta vacuum effects. 



4. Summary 

The QCD Lagrangian with 6 vacuum has the sign problem because of the topological term. 



This term can be vanished by the chiral transformation ( |1.5[ ), but the transformed Lagrangian 
has theta dependence in its light-quark mass terms. Using the fact that the P-odd mass has a 
much smaller scale than Aqcd, we have proposed a way of circumventing the sign problem. The 



reweighting method defined by ( |1.9[ ) may allow us to do LQCD calculations and get definite results 
on dynamics of 6 vacuum. 

Furthermore, we have investigated effects of the theta vacuum on the phase diagram for the 
realistic 2+1 flavor system, using the three-flavor PNJL and EPNJL models. Particularly in the 
EPNJL model that is more reliable than the PNJL model, the transition becomes first-order even at 
6 = when 6 is large. This result is important. If the chiral transition becomes first order at /I = 0, 
it will change the scenario of cosmological evolution. For example, the first-order transition allows 
us to think the inhomogeneous Big-Bang nucleosynthesis model or a new scenario of baryogenesis. 
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